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O . Abstract 

o 

^SJ ' Atmospheric diffusion of high energy cosmic rays is studied analytically and the obtained 

^ , integral electromagnetic fluxes are compared with the data measured by emulsion chamber 

O ' detectors at mountain altitudes. We find a good consistency between them when the average 

, nucleon inelasticity coefficient is varying between 0.50 and 0.65. 

m 

^ 1 Introduction 
> , 

^ The inelasticity is a global parameter of hadronic interactions and it is a crucial attribute 

t:J" ■ in the calculation of the fluxes in cosmic ray cascades. It is defined as the fraction of energy 

\ giving up by the leading hadron in a collision induced by an incident hadron on a target nucleon 
or nucleus. The inelasticity, Km, in a nucleon- induced interaction is related to the spectrum- 
weighted moments evaluated at 7 = 1 (where 7 is the power index of primary spectrum) as 



(N 
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1- Zmn{i = 1) = / V fNN{ri) drj 



p> il Several authors have suggested that the average inelasticity-coefficient is an increasing func- 

1^ ■ tion of the energy [1, 2], whereas others proposed that it is a decreasing function [3, 4, 5]. The 

behaviour of high energy cosmic rays, which reflects the nuclear interactions in the energy region 
I ~ 100 TeV, indicates a constant value of the mean inelasticity equal to 0.50 [6]. A decreasing 
\ with energy of this parameter points out a strongly suppression of the pion production in this 

j_j ■ energy region and does not agree with cosmic ray observations. 

d I The objective of this paper is to estimate the average inelasticity coefficient through the 

comparison of our calculations with the integral electromagnetic fluxes measured with large 
Emulsion Chamber Experiments at mountain altitudes. The main purpose is to describe in 
qualitative detail single events, i.e., shower cores in their early stage of development observed 
at mountain altitudes. 

Our calculation is performed by solving analitically the cosmic ray diffusion equations for 
electromagnetic showers in earth's atmosphere. 

In these calculations we adopt to the particle production an approximated Feynman scaling 
in the fragmentation region. We have used the data on proton targets and on nuclear targets 
[7] to obtain the pion spectrum weighted moments. Although exist several functional forms to 
fit the behaviour of rising cross-section we assume the following one 

a = ao{E/er ■ 

So, beginning with the measured all-particle primary spectrum we diffuse the particle showers 
down to the detection depths and we investigate our numerical results as a function of the 
adopted average inelasticity coefficient. 
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2 Hadronic and electromagnetic cascades in the Earth's 
atmosphere 

The differential hadronic fluxes, H{t, E) at a depth t in the energy range E and E + dE are 
obtained integrating the nucleon and pion diffusion equations using the semigroup theory. We 
assume that only the pions are generated in the multiparticle production neglecting the particles 
with small fractions such as kaons, heavy mesons, etc.. In this case the solutions are [8]; 



°° / lyj / + \nn 



(1) 



with 



Zmn{i, a) = [' h^iv) dv (2) 

Jo 

where 77 is the elasticity coefficient of the nucleons in the atmosphere that is distibuted according 

to /jvAr(?7). 

The pionic fluxes are 
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Z„^{j,a,m,n) = f .U^{x) dx (4) 
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The fN-n{x) and fmri^^) are the energy spectra of the pions produced in the nucleon-air and 
in the pion-air interactions respectively, and x is the Feynman variable {x = 2pl/VS). 
The solutions (1) and (3) are subjected to the boundary conditions 

N{0,E) = No E-^-^+^^ (6) 

and 

n(o,E) = o (7) 

that represents the nucleon and pion fluxes at the top of atmosphere. 

The above solutions are obtained considering the interaction mean-free paths in the power 
form, A niE) — Aq^ E^" {h = n or tt). 

Using naive approximations [9] we can write the solutions (1) and (3) in the following forms 

N{t, E) = N{0, E) e-*'^N{E) ^g-, 

and 

U{t, E) = N{0, E) ^\ , f e-*"'(^) - e-*"^(^)) . (9) 

u)n\E) - LO-„{E) V / 

The Wn{E), u)Tr{E) and the (jj/v7r(7, E) functions are described in the appendix A. 
The production rate of 7-rays originated from the 7r° decay is given by [10] 

roc fjj^f , . ^ 

P,iz,E,) = 2j^ _(p--- + p;^-ai'-) (10) 
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where and P^j. are the production rate of the charged pions from the nucleon-air 

and pion-air interactions and are represented by 

TjH-air H{z,Eo) , . dx 

In the expression (10) we have used that the secondary pions are formed with equal multi- 
pUcity, so the number of tt" is half of the number of charged pions. 

The differential flux of the electromagnetic component can be written as 

/>oo 

F^^e{t,E)= dz dE^ P^{z,E^) {e + ^){t- z,E^,E) (12) 
Jo Je 

where {e + j){t — z, E^, E) are the photons, electrons and positrons produced by a single photon. 
They are computed using the approximation A [11] and are written in the following form 

(e + 7)(t - z, E^, E) = ^J (^^) ^ (Ari(s)e^i + iV2(s)e^= W^*"^)/^") . (13) 

The functions -/Vj(s), Xi{s) and Xq are the parameters with standardized definitions in cas- 
cade theory [11]. These functions and the evaluation of the electromagnetic fluxes are presented 
in the appendix B. 

In order to compare our calculations for electromagnetic showers with the emulsion chamber 
data measured at mountain altitudes we need to obtain the integral fluxes of this component. 
They are, 

F^,e{t,>E)= / F^,e{t,E)dE (14) 
Je 



3 Comparison with experimental data 

Emulsion Chamber Experiments at mountain altitudes arc very simple detectors constituted 
of multiple sandwich of material layers (Pb, Fe, etc.) and photo-sensitive layers (X-ray films, 
nuclear emulsions plates, etc.). They can detect electromagnetic showers of high energy {Ee^-y > 
1 TcV) and due to its powerful spatial resolution and energy determination they can study in 
details the cosmic-ray interactions. 

The experimental data are taken from J. R. Ren et al. [19] for Mt. Kanbala; C. M. G. Lattes 
et al. [13] for Mt. Chacaltaya and M. Amenomori et al. [21] for Mt. Fuji. 

The experimental fluxes at Mt. Chacaltaya are revised according to the energy determination 
introduced by M. Okamoto and T. Shibata [22] where the LPM and spacing effects were taken 
into account. This revision followed the same procedure made by N.M. Amato and N. Arata 
[20] for determination of total hadron fluxes. 

In order to make a comparison with electromagnetic fluxes measured with these detectors it 
is necessary to check the main assumptions made in our calculations; 
a) Primary cosmic-ray radiation 

The energy spectra obtained by Emulsion Chamber Experiments at mountain altitude are 
in the region of 3 ^ 50 TeV for (e, 7) component and 3 ^ 40 TeV for hadronic shower. So, we 
need to include in our calculations the primary spectrum up to knee (« 10^'"^ eV). In the present 
work we assume the primary cosmic-ray radiation in integral form reported by J. Bellandi Filho 
et al. [14] from JACEE data [15]. 

iVi(0, >E) = 1.12 • 10"^'(i;/5TcV)"^ ' (/cm^.s.st) (15) 

where £'(TeV/nucleus). 

Nucleonic intensity at f = 0, equivalent to primary cosmic-ray radiation is given by 

N{{),>E)=A^-'^Ni{Q,>E) (16) 
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where A is the average mass number of cosmic ray nuclei; we adopt A = 14.0 which is the value 
in low energy region {E/pait < 10^^ eV). 

b) Collision mean free paths 

We use to fit the behaviour of rising with energy of the p-air inelastic cross-section the 
following form [7, 14, 16, 17], 

aNiE)=aNo{E/TeYr (17) 

with CTjVo = 300 mb and a = 0.06 [7, 14]. 

So, the collision mean free path of the nucleon in the earth's atmosphere is expressed by 

XN{E) = XNo{E/TeV)-" (18) 

with Ajvo = 80 g/cm2 [7, 14]. 

For the pion mean free paths we assume that Xt^^/Xno ~ 1.4 and that they have the same 
energy dependence, like in the nucleon case [14]. Thus, we have 

X^{E) = X^„{E/TeY)-" (19) 

with = 112 g/cm^. 

c) Secondary pions: 

The functions fNTt{x) and f-KTrix) are the energy distribution of secondary nucleons and pions 
for nucleon-air and pion-air interactions, respectively. We have used in this calculation the data 
on proton targets and on nuclear targets [7] to obtain these functions and they are written as 

U^{x) = 2.08 e-^^ (20) 

X 

these functions are for two charged states of pions and the contribution of leading pions are 
taken into account in the last one. 

Figures 1-3 show the comparison of our solution with the integral electromagnetic fluxes 
measured at Mt. Kanbala (520 g/cm^), Mt. Chacaltaya (540 g/cm^) and Mt. Fuji (650 
g/cm^), respectively. Two curves are also drawn in the figures for < Kn > = 0.50 and 0.65, 
taking into account the above mentionated itens. 



4 Conclusions 

We have solved the diffusion equations of cosmic-ray hadrons analytically and after this we 
derived the integral electromagnetic fluxes using the approximation A of the electromagnetic 
cascade theory. 

These fluxes at mountain atmospheric depths decrease when we include in our calculation 
the rising of the cross-section and the decreasing of the mean nucleon elasticity. 

We conclude that the naive assumptions for the secondary production based on approxi- 
mate scaling in the fragmentation region is possible to explain experimental data on integral 
electromagnetic fluxes detected in large emulsion chamber experiments. 

This scenario requires a mean nucleon inelasticity coefflcient varying from 0.50 to 0.65. 
At low energy (few TcV) the agreement between experimental data and the calculated fluxes 
shows that < Kn >~ 0.50, that is the adjust indicated by former workers in the cosmic ray 
literature. For E > 10 TeV we see from the figure that < Kj^ >~ 0.65. Therefore, with this 
simple analysis we see that the nucleon inelasticity coefficient appears to show a tendency of 
rising with energy according to Mini-jet Models [17] and Quark-Gluon String Models [2]. These 
integral fluxes are somewhat affected when the coefficient a changes in the interval 0.06 to 
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0.10, which is the appropriated values to explain the experimental data for mountain emulsion 
chambers experiments. In the latter the fit between our calculation and the experimental data 
is obtained for < Kn > changing between approximately 0.45 and 0.60. 

These results are in agreement with that of L. Jones [18] based in an analysis on inclusive 
reactions data from accelerator, with analytical calculations on nucleon fluxes at sea level [9] 
and with hadron and electromagnetic fluxes at mountain altitudes [8, 10]. 

Appendix A 

Hadronic cascades 

The functions com{E), u}.^{E) and Zj^^{'y,E) in the hadronic fluxes (8) and (9) are [9] 

1 / 1 „ , . 1 



and 



where 



and 



Wr,{E) = wr, = j-^ (^1 - 2^~^(7) - ^Z^N^, a) j (22) 

u^iE) = uj. E° °^ = (^1 + iz,,(7) - lz^A^,a)j (23) 

u^ME) = = 1^ InA^) dx (24) 

z^^{^,a)^ j n^-'^ fMNin) dn (25) 

Jo 

Z^^i^,a)= [ x^-"f^^{x)dx . (26) 
Jo 

The functions fN-wix) and f-Kvix) are the energy distribution of secondary nucleons and pions 
for nucleon-air and pion-air interactions, respectively (see item c of section 3). 
The nucleon elasticity distribution is assumed to be [9] 

hN{v) = {l+P)il-vf (27) 

where r] = E/E' and /3 fulfills a consistency relation between the mean elasticity and average 
inelasticity, so that <7V> + < K > =1, for nucleons. When /? = we have the 
nucleon uniform distribution of elasticity (< rj > = 0.5). 

Appendix B 

Electromagnetic cascades 

The production rate of 7-rays produced from 7r° decay expression (10) can be written as 



Pj{z,E^) = 



'y+l-a + zawN \ Ajvg A^o 



uin{E) - uj^{E) J +1- a + z au^) {uJi^{E) - uj-^{E)) A^o 

The Z-factors and Z^^,^ is the last expression are 



ZnAi - a + ^ a w^v) = / fN^{x)da 
Jo 



(28) 



(29) 
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and ^ 

Z^^{'^-a + zauJi)= [ x'<-°'+^°"^' f^^{x) dx i = n oitt . (30) 
Jo 

The expression (28) was obtained considering the expansion 
1 1 + a ln{E/Eo) 



^Oh 



H = N 01 TT (31) 



The differential flux of the electromagnetic component (12) is obtained from this production 
rate and evaluated at the poles si = j — a + z a iVf^ and S2 = j — a + z a co,^. 
The integral flux is put in the form 

F^,e(i, >E)= F^,,^{t, >E)+ F^,^e{t, > E) (32) 

with 

F^,,,{t,>E) = / dz (iVi(si)e^i(«i)(*-^)/^'' +7V2(si)e^=(^^)(*-^)/^«) Ci(si) 

Jo Si[Si + 1) V / 

(33) 

and 

F^,^,{t,>E) = / dz U(52)e^^(^=)(*-^)/^° +iV2(s2)e^^(^^)(*-^)/^°) ^2(^2) 

Jo S2(S2 + 1) V y 

(34) 



where 



and 



with 



and 



(a;.(£)-a;.(£;)) Ao, 



-^Ar,r(Sj) = / a;"' fNn{x) dx (37) 



x^' /,^(x) dx . (38) 



■^7r-7r(5i) — / 

Jo 



The functions Ni and Ai are the usual parameters with definitions in cascade theory 
and can be obtained in the reference [11]. Xq is the radiation length in the atmosphere 
{Xo = 37.1g/cm2). 
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Figure Captions 



Figure 1 - Integral electromagnetic spectrum at Mt. Kanbala (520 g/cvo?). The full lines 

represent the calculated flux for < K]^ >= 0.65 and the dashed lines for < Kj^ >= 0.50. All 
lines are for a = 0.06. The experimental data are taken from J. R. Ren et al. [19]. 



Figure 2 - Integral electromagnetic spectrum at Mt. Chacaltaya (540 g/cm^). The full lines 
represent the calculated flux for < Kn >— 0.65 and the dashed lines for < Kn >= 0.50. All 
lines are for a = 0.06. The experimental data are taken from C. M. G. Lattes et aL(black dots) 
[13]. 



Figure 3 - Integral electromagnetic spectrum at Mt. Fuji (650 g/cm^). The full lines represent 
the calculated flux for<iirjv>=0.65 and the dashed lines for < Kf^ >=0.50. All lines are for 
a = 0.06. The experimental data are taken from M Amenomori et al. [21]. 
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